scalar DiNum =  pyrolysis->solidRegionDiffNo();
